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(57) Abstract 

It is described a discrimination 
procedure of a wanted signal from a 
plurality of cochannel interferents re- 
ceived by array antennas of GSM or 
DCS base transceiver stations. The 
procedure includes a phase for the es- 
timate of the number and arrival di- 
rections of the interferents, and of the 
wanted signal, followed by a spatial fil- 
tering phase in which the signals trans- 
duced by the sensors of the array are 
linearly combined among them through 
multiplication coefficients, or weights, 
organized in a vector w satisfying the 
two following conditions: A) Spatial 
filtering constrains the gain of the array 
in the ratio between wanted signal and 
noise, compared to the traditional use 
of a single sensor, so that the gain is not 
less than a properly selected threshold; 
B) it minimizes the ratio between the 
sum of interferents' powers and wanted 
signal power. 
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DISCRIMINATION PROCEDURE OF A WANTED SIGNAL FROM A PLURALITY OF 
COCHANNEL INTERFERING SIGNALS AND RECEIVER USING THIS PROCEDURE 



5 

Field of the invention 

The present invention relates to the field of the so-called "intelligent" array antennas, 
and more in particular to a discrimination procedure of a wanted signal from a plurality 
of cochannel interferents received by array antennas of base transceiver stations for 
10 cellular telecommunication and relative receiver. 
Background art 

The use in mobile radio environment of antenna consisting of one or more 
arrays of electromagnetic field sensors (array) is known, for instance, from the 
European patent application published under No. 0593822A1. This document in fact, 

15 describes, mentioning the first claim, "An arrangement of antennas for a base 
transceiver station including a plurality of antenna arrays, each array able to form a 
multiplicity of narrow radiation lobes, separated and partially overlapped in the 
azimuth plane, the arrays being positioned in such a way that all the lobes generated 
by the same give a substantial omnidirectional coverage on said plane; lobes 

20 formation means (beamforming) in azimuth and angular elevation for each said array; 
a plurality of radiofrequency transceivers each one to transmit and receive the signals 
relative to one or more channels; a switching matrix to connect each transceiver with 
one or the other array through said beamforming means; etc.". 

More in particular, a beam former, for example a Butler matrix, is associated to 

25 an antenna array. A beam former consists of phase shifting and adder circuits, used 
in reciprocal way both in reception and in transmission. Considering for instance 
reception, the beam former has a plurality of input ports for the signals coming from 
the sensors of the array , and a multiplicity of output ports, each one relative to a pre- 
set azimuth direction and corresponding to a particular combination of module and 

30 phase of input signals. A dual behaviour applies to transmission, where the same port 
selected during reception is used to transmit towards the mobile . The antenna array, 
and relative beam former, are therefore an essential part of a system capable to 
identify the direction of the signal transmitted uplink by a mobile, that the system 
follows condensing a narrow radiation lobe in which the power of the signal 

35 transmitted down-link towards the mobile itself is concentrated. This equals, to an 
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intelligent behaviour of the antenna, which deviates from the traditional utilization of 
antennas in the same sector of the technique. Thanks to the intelligent behaviour, the 
interference from cochannel channel is reduced and the re-utilization of the same 
frequencies in adjacent cells is made possible. The deriving advantage is 
considerable and consists in the possibility to increase the dimension of cells, at equal 
transmitted power, in low traffic areas or reduce the frequency re-utilization distance, 
increasing the number of carriers per cell in high traffic areas. 

The estimate of wanted signal arrival direction avails of means for the 
measurement of signals present at output ports of the beam formers and of a 
processor evaluating the above mentioned measures and selecting the best direction, 
for instance that one for which the signal level is higher. 

The antenna arrangement described in the mentioned application, overcomes 
the limits introduced by sectorial antennas in corner-excited cells, still widely employed 
in the mobile radio field, due to need of infrasectorial handovers as the position of the 
mobile part around common antennas varies, though remaining close to the same. 
The intelligent antenna system enables in fact a transceiver to be associated to any 
possible narrow radiation lobes, sunburst arranged on the azimuth plane, therefore 
the cell appears to the network as if it were equipped of omnidirectional antenna, but 
without the relative known drawbacks. 
20 Background art 

Notwithstanding all the features mentioned above, the intelligent antenna 
arrangement described in the mentioned application, fails to exploit the capability of 
an antenna array to shape the overall radiation diagram so that nulls are steered in 
the interferents' directions. For example, if interfering and wanted signal arrival 
25 directions are different but contained in the same beam, no discrimination can be 
made between them. This stands also if beam is selected, after demodulation, on the 
basis of wanted signal quality. 

In order to exploit the interference nulling capabilities of the antenna array a 
numerical (and not analog) beamforming has to be performed on digitalized signals 
30 from the sensors of the array. 

Numerical beamforming is nothing more then a linear combination of base band 
signals from the array elements using a set of complex coefficients w. 
The procedures for numerical beamforming mainly differ in the choice of the 
coefficient's set w and usually 2 different methods can be used (see [l] for a more 
35 detailed description and for a comparison based on system capacity): 
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1) AH cochannel user arrival directions are estimated and coefficients w are chosen 
with the constraint that the overall array gain is "nulled* in the interferent arrival 
directions. See also [2] for an overview of known art about beamforming and "null 
steering" 

5 2) Coefficients w are calculated so that the difference between the combined signal 
and a reference part of the signal is minimized. 
Both these methods for coefficients w calculation suffer from the fact that they 
optimize uplink receiver performance but the coefficients w are not always reliable for 
downlink transmission. 

10 For example, when uplink noise level is negligible and the interferent and wanted 
signal arrival directions are close each other, both methods try to reduce interference 
as much as possible. The result is often a very low level of the array gain in the 
direction of the wanted signal and thus only a small fraction of transmitted energy will 
be directed towards the user. 

15 Summary of the invention 

The scope of the present invention is to overcome the above mentioned 
drawbacks and to indicate a discrimination procedure of a wanted signal from a 
plurality of cochannel interferents received by array antennas of base transceiver 
stations for cellular telecommunication. 

20 To attain these objects, scope of the present invention is a discrimination 

procedure of a wanted signal from a plurality of cochannel interferents received by 
array antennas of frequency division (FDMA), or time division (TDM A), or mixed 
FDMA/TDMA multiple access telecommunication systems, re-employing a same 
frequency group in adjacent territorial areas, including an estimate phase of the arrival 

25 directions of the said interferents and of said wanted signal, and a successive phase 
of spatial filtering in which signals transduced by the relative sensors of a said array 
are linearly combined among them through multiplication coefficients, or weights, 
obtaining a reception signal cleaned from interferents; said phases being repeated for 
each one of the time slots, in the TDMA frame, in which the wanted user transmits, 

30 characterized in that said weights satisfy the two following conditions: 

A) the gain of said array in the ratio between said wanted signal and the noise after 
spatial filtering, compared to the traditional use of a single sensor, is constrained so 
that it is is not less than a properly selected threshold; 

B) the ratio between the sum of interferents' powers and wanted signal power is 
35 minimized, as described in claim 1. 
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The subject procedure can find useful application in the realization of a receiver 
for base transceiver stations of cellular telephone systems of the GSM 900 MHz, or 
DCS 1800 MHz type 

The occurrence of condition A) enables to obtain some advantages: 
5 1. to improve the performance in uplink reducing the performance sensitivity of the 
spatial filtering to the error made in the estimate of wanted signal arrival directions. 
2. To enable, in downlink, a high directivity of the antenna array radiation diagram in 
the direction of the wanted. 

The occurrence of condition B) enables to improve the performance in the two 
10 connections, uplink and downlink, reducing the interference. 

Profitably, it is possible to express in the vectorial form the equations 
concerning the estimate phases of the arrival direction and spatial filtering, as well as 
those generated by the imposition of the above mentioned featuring conditions, thus 
making the mathematical computation for the determination of said weights more 
15 compact. 

Further object of the present invention is a frequency division multiple access 
receiver for telecommunication systems (FDMA), or time division multiple access 
(TDMA), or mixed (FDMA/TDMA), which re-employ a same frequency group in 
adjacent territorial areas, characterized in that it includes means for the actuation of 

20 the discrimination procedure of a wanted signal from a plurality of cochannel 
interferents, which already formed the object of an invention, as described in claim 13. 
Brief description of drawings 

Additional scopes and advantages of the present invention will be more evident 
from the following detailed description of an embodiment of the same and attached 

25 drawings, in which: 

- Figures 1 and 2 show an array of sensors used in the receiver of the present 
invention, when it is invested by plane waves coming from different directions; 

- fig. 3 shows a block diagram summarising the operational phases of the 
discrimination procedure object of the present invention; 

30 - fig. 4 shows the diagram of a function enabling to obtain a fictitious parameter 
useful to the calculation of weights entering the spatial filtering made according to 
the procedure of the invention; 

- fig. 5 shows a curve of the merit parameter G n (6) (Array gain) versus the arrival 
direction on the azimuth plane, for a receiver employing a discrimination algorithm 

35 according to the known art; . 
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- fig. 5' shows a table comparing the main merit parameters of the receiver scope of 
the present invention, in which it is included G n (9), with the corresponding 
parameters of receivers employing different discrimination algorithms according to 
the known art; 

5 - figures 6 and 6' represent a curve and a table, respectively relative to the same 
merit parameters of figures 5 and 5' calculated, contrarily to the first ones, 
employing reduced rank matrices; 

- fig. 7 shows a curve of the merit parameter G n (9) versus the direction on the 
azimuth plane for the receiver employing the discrimination procedure object of the 

10 present invention; 

- fig. T shows a table making the comparison of the main merit parameters of the 
receiver scope of the present invention, with the corresponding parameters of 
receivers employing the discrimination algorithms according to the known art, said 
parameters being calculated using matrices differently reduced compared to those 

15 of figures 6 and 6'; and 

- fig. 8 shows a general block diagram of the receiver scope of the present invention. 
Detailed description 

Making reference to fig. 1, we notice an array antenna consisting of N 
electromagnetic field sensors identified a 1f a 2> a, a Nt of known type, arranged 

20 in straight line and separated one from the other by a distance d, typically d = A/2, 
where X is the wave length of the radiofrequency in the middle of the band used by 
the particular mobile radio system employing the procedure scope of the invention (16 
cm approx. in case of GSM 900 MHz, 8 cm in case of DCS 1800 MHz). The array is 
invested by a given number of plane waves s 1t S/, .... s Ri whose arrival directions 

25 on the azimuth plane 9i, .... 9/ are indicated for two of them, Values 9 1t .... 9/ 
corresponding to angles formed by radiuses s^ t .... S/ with the line of sensors ai, a N . 
In a real scenario where the array is used by a base transceiver station (BTS) of a 
cellular telephone system, for instance GSM or DCS, plane waves Si, .... s y 
correspond respectively to echoes of a wanted signal (identified by index / = u) 

30 transmitted on an assigned carrier and time slot, and to interfering echoes generated 
in neighbouring cells by communication channels which re-employ the same 
frequency in the same time slot. 

Fig. 2 is useful to evaluate the progressive delay, or the corresponding phase 
shifting, according to which a plane wave s f invests the different sensors of the array. 
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If s(t) is a signal impinging on a uniform linear array from direction 9 t the signal 
received by the n m sensors x n (t) is a phase shifted replica of the s(t) according to the 
relation: 

x rt (i) = s(t)e J x withn=1...N (1) 

Thus the replicas on different sensors can be described by 



<(<)= 



"*,(')" 




1 

2arf(n-l)eos(e) 




"<(e)" 




*.(') 




e J - 

2nrf(A/-l)eos(6) 

/ x 


5(/) = 




s(t) = d-(e)s(/) 








am 





where ^*(9)is the response of the /? th sensor to a signal coming from the direction 9. 

Defining the response d*(6) as the complex conjugate of d n (Q) is useful, in the 
following, to write linear combinations as scalar products. 

In case of more signals Sj(t) impinging from directions Q u the signals x n (t) received by 
different sensors are described by: 



'(<) = 



*»(') 



(2) 



where n n (t) are the uncorrelated noise components on array sensors 

Making reference to fig. 3, we notice four blocks in cascade called ARY, 

STDOA, STPESI, and FSPAZ, respectively. The ARY block is the sensor array of fig. 

1 and 2 which gives the signals x n (t) given by the expression (2). 

The STDOA block processes the signals x n (t) to estimate the number R of plane 

waves s f that invest the array ARY and the relative arrival directions (DOA Directions 

Of Arrival) on the azimuth plane, given by angles Q t . The operation of the STDOA 

block is known, for instance from [3] and [4] . 

The information obtained by the STDOA block, in particular the directions e it 

where implicitly a particular value / = u of index / indicates the direction 6u of the 

wanted signal, is transferred to the next STPESI block that employs it in the 

calculation of appropriate spatial filtering coefficients w n , according to a novel method 

which shall be shortly described. 

Weights w n , together with the signals x n (t) transduced by the array ARY, are allowed 
to reach the FSPAZ block that gives a reception signal r{t) cleaned from interferents, 
thanks to a spatial filtering (beam forming) employing the following expression: 
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where it can be noticed that the spatial filtering is a simple linear combination of 
signals x„(t), transduced by the array ARY, weighed by coefficients w„. 

The (3) can be written in a form suitable to highlight the s, signals impinging on 
5 the sensors of the array ARY: 

r(t) = Z*(ofz*v£(e,)| + Z w 

/ V n ' n 

whose usefulness is that to enable the calculation of weights w„ according to the 
procedure scope of the present invention. It must be pointed out that incident signals 
s,<t) become accessible only through the corresponding transduced ones x„(t).. 
10 For the determination of weights w n it is convenient to rewrite the (4) in vectorial 

form highlighting the contributions of the wanted signal s u (t), of interferents Sj(t), and 
of the noise n(t), as it results from the following expression: 

r(t) = s u (tfd« w) + Z (*» + wn « (5) 

i 

- where: simbol H indicates a transposed and conjugated matrix 
15 - w is the N element vector of coefficients w n ; 

- is the N element vector containing the response of the array in the direction 0 U 
of the wanted signal s u (f ) ; 

- cf is the N element vector containing the response of the array in the i th interfered 
direction ; the index / of the summation extends to value R-1;n(/) is the N 

20 elements noise vector. 

As for the operation of the STPESI block, it is convenient, with the aid of the (6), 
to introduce two gain parameters, respectively identified G n and G* in order to be able 
to introduce constraint conditions on said parameters for the calculation of weights w n 
. The above mentioned gains are the following: 

25 Gain of wanted signal on noise: 



G (6) 



llwll 2 



is the gain in the S/N ratio versus the traditional case of use of a single element 
antenna. The notation (...[indicates the Euclidean norm. The maximum value of G„is 
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equal to |d| 2 

Gain of the wanted signal on the interfered h 



o.-Fr=!r P> 



icfw 



is the gain in the C/l ratio, concerning an interferent /, versus the traditional case of 
using a single element antenna, Analogously the gain of interferent on wanted signal 
is given by: 



1 llcfwfi 2 



With these assumptions, the STPESI block calculates the vector of weights w 
satisfying the two following conditions: 
10 A) constrains the gain G„ in order that it is higher than, or equal to, a fraction a 2 of the 
maximum value, in order to obtain the above mentioned advantages. The 
parameter <x 2 can be varied in a continuous fashion in order to trade off 
interference reduction for directivity towards the user The condition A) is equivalent 
to the expression: 

B) minimizes the sum of gains of interferents on wanted signal, maintaining the 
constraint A) to the purpose of reducing the interference in uplink and downlink. 
Condition B) is equivalent to the following expression: 

f|Cwf 



mm 



n2\ 

(9) 



20 where C" is a matrix formed by vectors c" of the (5). 

Should all the interferents have the same power the (9) is equivalent to minimize the 
ratio between the sum of interferents' powers and the power of the wanted signal 

(//C). 

The STPESI block calculates the vector of weights w imposing the conditions 
25 (8) and (9). To simplify the computation and without involving any conceptual limitation 
to the procedure, it is imposed that: 
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d"w = l (10) 
The above mentioned normalisation condition involves the multiplication of vector w 
by a constant scale factor, possibly complex, not influencing the spatial filtering made 
by the equation (5) since, as it can be noticed from the same, the wanted signal, 
interferents and noise would be multiplied by the same constant, and therefore the 
filtered signal, which can be denormalized dividing it by ( the same scale factor. With 
the setting of the (10) the (8) becomes: 



llwll 

and the (9) becomes: 



in(||C*w|| 2 ) O') 



10 mini 

The STPESI block employs the method of Lagrange multipliers to make a 
constraind minimization of (9'). To this purpose, it is first constructed a function 
F(w,X,p,x) of variables w.X.p.x , indicated with (11): 

15 F(w. X, 0, X ) = |C* wf 4- ^«wf 2 - + 2p(Re(d w w) - 1) + 2x(lm(d" w) - l) 

where parameters w,X.,p,x are arbitrary constants to be determined equalling the 
gradient of F(w,X,p,x) to zero. The resulting equations are given by: 

3F(w.X,P,x) = 2CcH w + 2 ^ w + 2d / p + Jy \ = 0 

aF(wA,p,x) = |wf _i =0 (12) 

dX a |dj m 



r aF(wA,p >x ) .^(w.x,p, x )l 2(d „ w _ l) = 0 
I ap J d% ) 1 ' 



Solving for p.x and replacing it is possible to obtain the value of vector w according to 
the X parameter, the resulting expression is the following: 

w _ (cc*+xi)-'d 

d"(CC"+Xl)~'d 
25 where I is a unit matrix of the same order of the matrix C. 
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To solve the (13) it is first necessary to determine the value of the parameter X. 
This real scalar is obtained assigning the desired value of the gain G„ given by the 
(8*), that is equivalent to remove the > sign from the same expression, having 
imposed a precise value of a 2 . The (8*) therefore becomes: 

Using a known factorization of eigenvalues and eigenvectors of the CC" 
matrix of the type: 

CC" = E"AE ( 14) 
where A is the diagonal matrix of the eigenvalues of CC", and E is the matrix of 
eigenvectors, we obtain the following expression of the gain (8"): 



j Ml 

kf = T e; — (15 > 

where A. are the eigenvalues, real and not negative, of matrix CC" ; e s are the 
elements of the vector e = jE"d| more particularly, the elements e, correspond to the 
module of the components of vector E"d ; and X is the above mentioned real scalar, 
considered continuously variable to suit the arbitrary character of the selection of a-. 
From the (8") and (15) the following equation is reached in the unknown X : 



1 VT A + 



whose solution gives the numeric value of parameter X that replaced in the (14) gives 
20 in its turn the vector of weights w. 

In view of a reduction of the computational complexity requested to the STPESI 
block, once the numeric value of the unknown X is obtained as said above, these can 
be replaced in the following equation: 



WO 99/33141 



11 



PCT/EP98/0816S 



E(A + 7U)''E*d 
W d"E(A+-Al)- , E w d (17) 
obtained replacing the (14) in the (13). The (17). contrarily to the (13). has the 
advantage not to require difficult matrix inversions, except for the easy inversion of the 
term in brackets which is a diagonal matrix. 
5 As far as the solution of equation (16) in the unknown value X, fig. 4 shows with 

a continuous line the diagram of the function j^jrft) (or equivalency cr(X)) and 

with a dotted line the diagram of the function to be maximised ? (X) that 

|C*w| 

corresponds to the gain in C/l in case of equal power interferents. 

In fig. 4 the x-axis indicates the X value normalized at the value of the maximum 
10 eigenvalue max(Aj) while on the y-axis values are in dB . 

As it can be noticed from f,g.4 the function is always increasing for X>0 

and tends to 1 for X-*» 

The study of both functions suggests a procedure for the solution of the (16) in 
the unknown X . To this purpose the following steps are foreseen: 

" wecalcu,ate MW (x=0); 

~ ' f flwlPfld! 2 (A = 0) " a2th6n We Ch °° Se X= °' because « is already the desired value 
since for X = 0 we obtain the unconstrained minimum value of |c"w| 2 ; 



- if 



1 



Mffdi 2 ^ " °) K a * then X is increa sed from 0 employing an iterative method for 

equation solution. In particular using the first order Newton method we obtain the 
iterative relation: 
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where function f(x) is properly selected in order to guarantee the sequence 



convergence. Practically f(x) can be chosen as f(x) = log(x) or f{x) = 



x-1 



x + 0.05 

For example from the diagrams of fig. 4 it results that the case represented 
5 corresponds to a very low noise gain G n for X = 0. On the contrary, selecting a 2 = 0.5 
(-3dB) so that the gain in S/N is only 3 dB below the maximum achievable value, it 
results that the solution of (16) is approximately X = 0,15xmax(A). 

Referring to figures 5, 5', 6, 6\ 7 and T the performance of the discrimination 
procedure scope of the present invention are now compared with similar 
10 performances of the main beamforming algorithms known so as described in pj. The 
scope is that to highlight additional drawbacks of the known art now superseded by 
the present invention. 

Merit parameters, expressed in dB, highlighted in the tables of figures 5\ 6' and 
T attached to the relative figures 5, 6 and 7 are the following: 
15 - the gain G„ in the S/N ratio, versus a traditional antenna having a single element, 
given by (6); 

C 

- the gain in the — ratio between the wanted signal and the interferent signal at the 

output of the beamformer (FSPAZ block), versus a traditional antenna having a 
single element; 

20 - the gain in the — - ratio between the wanted signal and the noise plus 

the interferent signal at the output of the beamformer (FSPAZ block), versus a 
traditional antenna having a single element. Concerning this last parameter, the 
following expression applies: 
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Using the same terminology of [2] , the known algorithms considered are the 
following: 

1. Minimum Variance (MV): known the wanted signal direction, the weights w n are 
selected in order to minimize the power at the beamformer output with the 

5 condition d w = 1 . 

2. Linear Constrained Minimum Variance (LCMV): in addition to the condition 
d"w = 1 , the weights w n minimize the power at output of the beamformer with the 
additional constraints that the array response is null in the interferents directions 
cfw = 0 . 

10 3. Linear Constrained (LC): weights w n meet the previous constraints on the wanted 
signal and on interferents, also minimizing ||w|| 2 . 

The procedure scope of the present invention is indicated in the tables with: 
4. Constrained Gain Minimum Interference (CGMI): weights w n minimize the sum of 

gains of interferents on wanted signal (fc* wf j , with the constraint that the gain G n 

15 of the wanted signal on the noise is higher than a pre-set threshold (6 dB in the 
following examples). 

To obtain the above mentioned tables and comparison diagrams the following 
scenario was selected, particularly suitable to highlight the defects of the different 
algorithms: 
20 • an array of 8 elements is used 

• a wanted signal with arrival direction orthogonal to the array plane (DOA = 0°); 

• two interferents of equal power, one half of the wanted signal , with a small angular 
separation from this last: (C/l = 0 dB; DOA = [-5°,5°]); 

• very low noise level (S/N = 40 dB); 

25 • arrival directions of the wanted signal and interferents perfectly known. 

Making reference to the comparison table 1 of fig. 5', we can notice that the first 
3 algorithms show poor performance due to a really low gain G n (-25 dB for LCMV and 
MV and -5 dB for LC), contrarily to what takes place for the CGMI method, which has 
a constrained G n equal to +6 dB, used in the discrimination procedure scope of the 

30 present invention. This means that, when using the LCMV or MV algorithm, the power 
radiated in downlink towards the mobile is 25 dB lower compared to the case of single 
element antenna. In this condition, also if interference is absent, the S/N at the mobile 
could be too low for a correct detection. 
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In the corresponding fig.5 t showing the radiation diagram that is the gain g„(9) 
of the array concerning the application of the algorithm LCMV, we notice that the 
above mentioned diagram shows a very low value (-25 dB) in the wanted signal 
direction, null (-00 dB) in the directions of interferents, and high values in the other 
values of 6. Thus fig. 5 shows how, with algorithms from known art, most of the power 
is radiated in directions different from the wanted direction. 

This behaviour is due to the fact that the LCMV algorithm, as the other known 
algorithms can run up against the inversion of almost singular matrixes, that is with 
very high condition number, when arrival directions of interferent and wanted signal 
are close each other. 

A way to mitigate this effect is to use matrixes of reduced rank, obtained 
selecting only the highest singular values of the matrixes to invert. 

Figures 6 and 6\ maintaining the previous scenario, show the effects obtained 
employing reduced matrixes in the different algorithms indicated; the reduction is 
conducted selecting only the 2 highest singular values. Concerning the CGMI method 
used in the invention, it does not involve the inversion of nearly singular matrixes, 

The radiation diagram of fig.6, still showing the array radiation diagram g„(8) 
concerning the application of the algorithm LCMV, has a much more higher directivity 
in the wanted, however, also interferent signals result very high. Table 2 of fig.6' 
confirms this behaviour and shows that with the known algorithms a high g„ value is 
obtained in the wanted signal direction, but the capability to null the interferents is lost, 
as it results from the low values of parameter A(C/I). On the contrary, with the CGMI 
method, in which the gain G n is previously set at 6 dB, a higher reduction of 

interferents is obtained. 

Figures 7 and 7\ maintaining the previous scenario, show the effects obtained 
increasing to 3 the number of the highest singular values used in the inversion, in 
order to partially recover the capacity to nullify the interferents, though accepting to 
lose something in the gain g„ since we come mQre close to the situation of figures 5 
and 5\ Actually, the results of Table 3 of fig.7' show the correctness of the reasoning, 
but show also that the gain g„ of known algorithms is too low for a correct antenna 
operation, contrarily to what takes place for the CGMI method used in the invention. 

Fig.7 shows the radiation diagram obtained applying the CGMI method used in 
the present invention, confirming the results shown in the table. 
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What we can conclude from all comparisons made is that, differently form 
known algorithms, the algorithm which is the scope of the present invention allows to 
easily and continuously trade off antenna gain in the wanted signal direction for 
capacity to null the interferents. The conclusion is supported by the fact that the 
5 passage from the situation of Table 2 (fig.6') to that of Table 3 (fig.7'), describing 
known art results, does not allow continuously variable intermediate situations, since 
reduced matrices can be only obtained selecting an integer number of singular values. 

Fig. 8 gives a general representation of a receiver for base transceiver station of 
a cellular telecommunication system GSM, or DCS type, employing the discrimination 

10 procedure scope of the present invention. More particularly, the receiver of fig.8, and 
the relative associated transmitter, can find application in frequency division (FDMA) 
or time division (TDMA), or mixed FDMA/TDMA multiple access systems re-employing 
a same frequency group in adjacent territorial areas. 

Making reference to fig. 8, we can notice the array antenna ARY consisting of N 

15 elements a 1t .... aj, .... a N , connected to a block RIC-FDMA/TDMA from which digital 
demodulated signals CHm, come out at each time slot, reaching a process module 
PROC. In the symbol CHmj, the index m indicates the generic carrier m-th of the M 
separate carriers assigned to the receiver, while the index j indicates the j-th of the N 
replica of the reception signal relative to the nrvth carrier. A snapshot is thus supplied 

20 of the signals coming out from the above mentioned block at the current time slot. The 
process module PROC processes, according to the procedure described above, the N 
replica of each one of the M channels CHm and supplies M digital signals CHFm, 
spatially filtered, corresponding to the M channels simultaneously received. 

In the operation, the RIC-FDMA/TDMA block performs all the operations 

25 necessary to the reception of M channels of the FDMA/TDMA type from each one of 
the N sensors of the array, that is: 

• radiofrequency filtering for the suppression of the spurious out of the total reception 
band and subsequent frequency splitting; 

• conversion at intermediate frequency and filtering; 
30 • conversion analogue to digital; 

• best demodulation and filtering to obtain demultiplexed signals FDMA/TDMA, to be 
sent to the PROC block. 

The M channels CHF1, .... CHFm, .... CHFM, spatially filtered, coming out from 
this last block suffer the following additional processing inside the receiver of fig.8: 
35 • reconstruction in base band of the M original transmission bursts associated to the 
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M channels and positioning inside frames and multiframes; and 
• new coding of digital bursts in a format (PCM 30 channels) compatible with the 
protocol (LAPD) adopted on beams connecting the base transceiver stations to the 
relative station controller (BSC). 
5 The original character of this receiver is of course that of the PROC block, so 

that the RIC-FDMA/TDMA block can be considered known to the Field technician. 
More particularly, it is possible to adopt for the above mentioned block, an architecture 
based on single narrow band receivers, as for instance results being that of the 
disposition of antenna according to the known art mentioned above, or a disposition 
10 based on the use of broad band transceivers simultaneously processing several 
channels, like for instance the one described in the European patent application No. 
97830229.7 filed under the name of the same applicant. 

Profitably, the PROC block of fig. 8 can be implemented through a 
microprocessor for the mathematical processing of digital signals (DSP), or more 
15 adequately, through digital integrated circuits of the ASIC type (Application Specific 
Integrated Circuit). 

Therefore, while a particular embodiment of the present invention has been shown 
and described, it should be understood that the present invention is not limited thereto 
since other embodiments may be made by those skilled in the art without departing 
20 from the scope thereof. It is thus contemplated that the present invention 
encompasses any and all such embodiments covered by the following claims. 
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CLAIMS 

1. Discrimination procedure of wanted signal (s u (t)) from a plurality of 
cochannel interferents received by array antennas (ARY) of frequency division 
or time division or mixed multiple access telecommunication systems, re-employing a 

5 same group of frequencies in adjacent territorial areas, including an estimate phase 
(STDOA) of the number of said interferents and of the their arrival directions and of 
the arrival direction of said wanted signal 0 Ul and a subsequent spatial filtering phase 

(FSPAZ, STPESI) in which signals x n (t) transduced by relative sensors fa a N ) of a 

said array (ARY) are linearly combined among them through multiplication coefficients 

10 w n , or weights, obtaining a reception signal r(t) cleaned from interferents; said phases 
being repeated for each one of the time slots, in the TDMA frame, in which the wanted 
user transmits, characterized in that said weights meet the following conditions: 

A) the gain of said array in the ratio between said wanted signal and the noise 
components after spatial filtering, compared to the use of a said single sensor, is 

15 constrained in order to result higher than or equal to a properly selected threshold; 

B) the sum of gains of each said interferent on the wanted signal after spatial filtering, 
compared to the use of a said single sensor, is minimized . 

2. Procedure according to claim 1, characterized in that said spatial filtering is 
made employing the following expression valid for complex envelopes in base band: 

n 

where: 

- r(t) is a filtered signal whose characteristics reflect said discrimination of the wanted 
signal from said interferents; 

- Xn(t) are the said signals transduced by the relative said sensors (a lt .... a N ); and 
25 - w n are said multiplication coefficients, or weights. 

3. Procedure according to claim 2, characterized in that said expression used to 
perform the spatial filtering is equivalent to the following one, given in vectorial form: 

K<)= s 9 (t%d H w) + E*i('X c N) + wn(/) 

i 

where: 

30 - s u (t) is said wanted signal inciding on-said array (ARY); 

- Si(t) is an interferent Mh of said plurality inciding on said array (ARY); 

- w is the N elements vector of said weights w n , relative to said sensors (a, a N )\ 
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- d" is the A/-element vector that indicates the response of the array in the direction 
of the wanted signal; 

- c" is athe A/-element vector that indicates the response of the array in the 
interferent direction hth of said plurality; and 

- n(/) is a noise vector at N elements, each one associated to one said sensor. 

4. Procedure according to claim 3, characterized in that said condition A) is 
equivalent to the expression: 



id* 



w 



w 



and said condition B) is equivalent to the expression: 



mm 



' |C"w|| : 



where: 



- a 2 is said pre-set fraction of said maximum value Id|| 2 = max 



f|dM 



being 



||...|| the Euclidean norm; and 

- C H is a matrix formed by vectors c" . 

5. Procedure according to claim 4, characterized in that said spatial filtering is 
made imposing the following normalization constraint: 

d"w = l. 

6. Procedure according to claim 5, characterized in that said vector of weights w 
is calculated through the following expression: 

(CC" -hXI^d 
d^CC^+Jd^d 
where: 

- I is a unit matrix of the same order of said matrix C ; 

- X is a real scalar parameter obtained solving the expression related to said 
condition A), in presence of said normalization constraint and imposing a precise 
value to said pre-set fraction a 2 . 

7. Procedure according to claim 6, characterized in that said matrix CC H is 



w = 
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submitted to a factorization at the eigenvalues and eigenvectors of the type: 

CC" = E"AE 

where: 

- A is the diagonal matrix containing the eigenvalues of said matrix , and 
5 - E is the matrix of the eigenvectors. 

8. Procedure according to claim 7, characterized in that said vector of weights w 
is also calculated through the following expression: 

E(A + XI)-'E*d 
d w E(A + M)-'E"d ' 

9. Procedure according to claim 7, or 8, characterized in that said parameter X 
10 is obtained solving the following equation having one unknown value: 



fiw| 2 



= a* 



where: 

- A, are the eigenvalues, real and not negative, of the matrix CC" ; and 

- e, are the elements of a vector e = |e "d| . 

15 10. Procedure according to claim 9, characterized in that the above mentioned 

equation in said unknown X is solved in iterative way through the following steps : 

- i. is catenated ^(^o); 

" ' f Mfldfl 2 ^ = °^ We choose X = °- becau se it is already the desired 

value since from the per X = 0 we obtain the minimum value of [c"w| 2 
20 ~ ' f Npjdp ^ = °^ ' X ' S increased from 0 usin 9 the iterative relation: 
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^+1 = *n + 



/(« 2 )-/(, 




where function f(x) is properly selected in order to guarantee the sequence 
convergence 

Iterations are executed until it results — ^ — r =ct 2 



1 1 . Procedure according to any of the previous claims, characterized in that 
N = 8 is said number of sensors (a,, .... a w ), and a 2 = 0.5 is said pre-set fraction of 
said maximum value ||d|| 2 . 

12. Procedure according to any of the previous claims, characterized in that 
said telecommunication system is a cellular telephone system. 

13. Receiver for frequency division, or time division, or mixed multiple access 
telecommunication systems, which re-employ a same frequency group in adjacent 
territorial areas, applying the discrimination procedure object of the previous claims, 
including: 

- an array of sensors of electromagnetic field (ARY); 

- radiofrequency filtering means (RIC-FDMA/TDMA) of the signal coming from each 
sensor (a u .... a N ) for the suppression of the spurious out of the global reception 
band, having downstream splitting means of a plurality of earners; 

- means for conversion and filtering (RIC-FDMA/TDMA) at intermediate frequency; 
analogue-to-digital conversion means; demodulation and filtering means to obtain 
in base band demultiplexed signals FDMA/TDMA of relative communication 



- process means (PROC) of said signals in base band and of reconstruction of the 
information originally transmitted on said channels; 

characterized in that said processing means include also: 

- estimate means (STDOA) of the arrival direction 9 U of a wanted signal s u (t), and of 
the number and arrival directions 6jOf a plurality of cochannel interferents signals 
si(t); 

- spatial filtering means (FSPAZ) of said signals in base band, to obtain a reception 




channels; 
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signal r(t) cleaned from said interferents, said means linearly combining among 
them x„(t) signals coming from the relative said sensors (ai, .... a N ) of a said array 
(ARY), through multiplication coefficients w n , or weights; 

- calculation means (STPESI) of the gain of the wanted signal on the noise, spatially 
5 filtered, compared to the traditional case of utilization of a single element antenna; 

- calculation means (STPESI) of the gain of the wanted signal on the interferent /-th, 
spatially filtered, relating to an interferent /-th and for all the interferents, compared 
to a traditional case of use of an antenna to a single element; 

- means calculating said weights (STPESI) imposing a value of said gain of the 
10 wanted signal on noise, higher than or equal to a duly selected threshold, and 

simultaneously a minimum value of the sum of gains of each said interferent on the 
wanted signal. 

14. Receiver according to claim 13, characterized in that N = 8 is said number of 
sensors (a^ a*), and a 2 = 0.5 is said predetermined fraction of said maximum 

15 value. 

15. Receiver according to claim 13, or 14, characterized in that it is used in a 
cellular telecommunication system. 

16. Receiver according to any claim 13 to 15, characterized in that said process 
means (PROC) consist of digital integrated circuits of the ASIC type {Application 

20 Specific Integrated Circuit). 
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